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The electron transport of different conical valleys is investigated in graphene with extended line- 
defects. Intriguingly, the electron with a definite incident angle can be completely modulated into 
one conical valley by a resonator which consists of several paralleling line-defects. The related 
incident angle can be controlled easily by tuning the parameters of the resonator. Therefore, a 
controllable 100% valley polarization, as well as the detection of the valley polarization, can be 
realized conveniently by tuning the number of line-defects and the distance between two nearest 
neighbouring line-defects. This fascinating finding opens a way to realize the valley polarization by 
line-defects. With the advancement of experimental technologies, this resonator is promising to be 
realized and thus plays a key role in graphene valleytronics. 

PACS numbers: 73.22.Pr, 73.61.Wp, 73.63.Bd, 85.75.-d 



I. INTRODUCTION 

Since the successful fabricating of graphene in 2004, 1 
its excellent transport properties) 2 ^ such as the excep- 
tional electron and thermal transport properties at room 
temperature, have attracted a lot of research interest to 
apply it in making future electronic devices EDEU Many 
ideas and attempts have been going on. One of the most 
fundamental and revolutionary attempts is to control the 
valley degree of freedom to realize the signal's storage and 
transport P^HIS] Supposed an achievement of this idea in 
future, graphene will probably replace all semiconductors 
in conventional electronics and bring a new era in many 
fields. Before the coming of the new era, there is a long 
way to improve the control over transport in graphene. 

Inspired by the crucial role of the doping in semicon- 
ductor, many researche rs ha ve studied the effects related 
to disorder in graphene lUilHUl w hich consists of adsorbed 
atoms (or molecules), charged impurities, vacancies, or 
other topological defects. It is anticipated to modify the 
transport property of pristine graphene to fulfill the re- 
quirement. Recently, a peculiar topological line-defect in 
graphene was reported experimentally by Lahiri et alPSl 
This topological line-defect is created by alternating the 
Stone- Thrower- Wales defect and divacancies, leading to 
a pattern of repeating paired pentagons and octagons, as 
shown in Fig. 1. In particular, it is found in this exper- 
iment that this line-defect has metallic characteristics. 
Great attention has been attracted to the line-defect and 
further investigations manifest a variety of its promising 
applications! 21 " 26 ^ 

Although a lot of theoretic proposals are presented in 
graphene to realize the valley filter before, all of them 
are not realized till now due to the extreme difficulty in 
the experimental implementation. After the discover of 
line-defect, Gunlycke and White^ reported that a val- 
ley polarization near 100% can be achieved by scatter- 



ing off a line-defect! 2 ^ This original idea subsequently 
arouse physics community's intensive interests! 27 " 29 * As a 
new discovery, the extended line-defect is very promising 
to realize the valley polarization in graphene. However, 
because the 100% valley polarization only maintains at 
a large incident angle in these reported results pMH] a 
serious problem will be encountered in experiment that 
the electron must always follow the direction of the line- 
defect to maintain a high valley polarization. This trans- 
port characteristic is almost fatal to utilize efficiently the 
high valley polarization. The purpose of this paper is to 
address a feasible approach to realize a 100% valley po- 
larization at a small incident angle by using a resonator, 
which consists of several paralleling line-defects. 

Given a high valley polarization by scattering off one 
line-defect, a question is putted forward naturally: What 
will happen in the presence of two or more line-defects? 
Motivated by finding an efficient realization of the high 
valley polarization, we therefore study the electron trans- 
mission coefficient in the presence of several line-defects. 
Intriguingly, it is found that an obvious change happens 
to the transmission coefficients with varying the distance 
between two line-defects. In particular, for the electron 
with a definite momentum and incident angle, the trans- 
mission coefficients can reach to one at a definite dis- 
tance between two line-defects. Simple numerical analy- 
sis argues that the large transmission coefficients should 
be attributed to the standing wave structure due to the 
electron's quantum interference, which is also named as 
Fabry-Prot interference in optics. Based on a further nu- 
merical simulation, an efficient device to realize a 100% 
valley polarization is further proposed for six line-defect 
configuration at the end. This fascinating finding opens 
a way to realize the valley polarization by line-defects. 

The rest of this paper is organized as follows: In Sec.pH 
we introduce the Hamiltonian in the tight-binding model 
and present the formula of the transmission coefficient. 
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Figure 1: (Color online) (a) The schematic diagram of the 
graphene sheet with a line-defect along y axis, (b) The Bril- 
louin zones of Graphene, which correspond the single cell and 
the supercell shown in (a), (c) The simplified lattice model 
of the infinite graphene with a line-defect. Here 
from the Fourier transformation along y axis. 
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Figure 2: (Color online) (a) The schematic diagram of K or 
K' Dirac cone . (b) The relationship between the electron's 
incident angle and momentums. (c) and (d) manifest the 
transmission coefficient in K and K' valleys. In this case, 
there is only one line-defect. Here, the dotted line represents 
the approximate result in the low energy limit: (1 ± sina)/2 



The numerical results are discussed in Sec. |III| Finally, a 



brief summary is given in Sec. IV 



II. THEORETICAL MODELS 

In the tight-binding approximation, a single layer of 
graphene with a line-defect, sketched in Fig. 1(a), can 
be described by the following Hamiltonian: 



H 



<iJ> 



where c\ and c iy ^ a /p represent the electron annihilation 
operators on the graphene site i and line-defect sites, re- 
spectively, and (•) refers to the nearest neighboring sites. 
Supposed a va riation less than 5% in the hopping terms 
T\ and T2 JHEmU j iere we a( j pt t\ « r 2 ~ t = — 1. 

Because of the existence of line-defect, the translation 
symmetry along x axis breaks down for a graphene sys- 
tem, as sketched in Fig. 1(a). However, the translation 
symmetry along y axis still retains and therefore k y is 
a good quantum number, namely a conserved quantity. 
For the sake of the convenience in computing, we do the 
Fourier transform only along y axis as follows: 



2-2 Jx>u 



% 2 ku 



c + 

i y ,a 



Y] C k y ,a 



i2kyi y a 



k v 



Cfc„ 



(2) 



a represents the lattice constant in pristine graphene. Af- 
ter the Fourier transformation, the Hamiltonian in Eq. 
(|TJ) is decoupled into H = H^ y . Note that here H^ y 
can be effectively represented by a quasi-one-dimensional 
lattice model, as shown in Fig. 1(c). Readily, we can de- 
scribe the Hamiltonian as follows: 

H K = ~ Yj ^1,1^1^2 - Yl ^1,2^2^,3 - Y vUTlVi,* 
iii 

~ Y ^1,4^2^+5,1 - (pl 1A f 2 (p ~ ^0^2^1,1 

-iplf 3 (p + h.c, (3) 

1 1\ ^ f 1 0\ ^ f 1 

1 



Ti = 



i20 i | j T 2 - ( q 1 1,^3 — 



where = k y a, <p[ a = [4 y ,*,a,A>4 v ^a,ij]> and f re P re " 
sents the unit length between the neighboring supercells 
at the graphene part. Here, i represent the position of 
a supercell, a takes the integer number from 1 to 4 to 
denote the different columns in a supercell, and A/B in 
4 i a a/b corresponds to the up/down site in the same 
column in Fig. 1(c). 

According to the Landauer-Biittiker formula and the 
Hamiltonian in Eq. (J3|, the transmission coefficient of 
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Figure 3: (Color online) (a) The schematic diagram of the 
graphene sheet with two line-defects. The transmission co- 
efficients in K (b) and K' (c) are mapped as a function of 
the incident angle a in the presence of two line-defects. Here, 
electron's energy is set as E — 0.01. 
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Figure 4: (Color online) Contour plots of the transmission 
coefficients in K (left) and K' (right) valleys. In this case, 
two line-defects are considered and the fixed distance between 
them is N d = 74. 



Besides, some previous papers reported that in the 
vicinity of the Dirac point, the transmission coefficient 
of low energy electron will obey the law: T K j K t = 
(1 ± siim)/2. In order to make this work essentially self- 
contained, we also confirm this conclusion and present 
the related derivations in Appendix. 



the line-defect system at zero temperature and low bias 
voltage can be represented asP^^l 

T ky (E)=Tr[Re(T L G r T R G a )} 1 (4) 

where G r l a is the retarded/advanced Green's function 
related to the line-defect Hamiltonian and T L / R = 
i^L/R ~ ^l/r) w ^ n t ne retarded/advanced self-energy 
Note that here the left /right lead is represented 
by a semi-infinite quasi-one-dimensional graphene lattice 
sketched in Fig. 1(c) and the sample Hamiltonian refers 
to the line-defect Hamiltonian, namely T3. 

It is worth noting that we can readily compute the 
transmission coefficients of the electron with a fixed in- 
cident angle a in two Dirac points, K and K' by using 
the above formula. As shown in Fig. 1(b), two Dirac 
points K and K' which are located at [0,±47r/3a] for 
a pristine graphene are now cast at ka and k' d for a 
graphene with a line-defect, [0, d=7r/3a]. Therefore, an 
electron in two valleys satisfies the following relations: 
k x — Qx and k y = ±7r/3a + q y: where q x /q y represents 
electron's group velocity along the x/y direction. Com- 
bined with the Dirac electron's linear dispersion relation 
£ _ Vsaq _ ^ a \J<z x + %j_ ^ transmission coefficients 
can be readily mapped as a function of the electron in- 
cident angle a where a = arctan^/g^). For simplic- 
ity in discussing, we will never distinguish the difference 
of the Dirac points in two Brillouin zones in the follow- 
ing and uniformly adopt K and K' to represent the two 
nonequivalent Dirac points. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. One line-defect 

Even though the transmission coefficient in the vicinity 
of the Dirac point was g iven by an approximated method 
in many references j 26 l 28 l to the best of our knowledge, the 
exact numerical calculation without any analytical ap- 
proximations is not unambiguously shown before. Given 
the realization of the valley polarization using the line- 
defect in future, what the transmission coefficient behav- 
iors at a relatively large energy is a fundamental question. 

In Fig. 2, the transmission coefficient is mapped as 
a function of the electron's incident angle for different 
energies. It can be observed the line of E = 0.001 (red) 
overlaps exactly with the dotted line, (l±sina)/2, which 
represents the transmission coefficient of low energy elec- 
tron in the vicinity of the Dirac point by an approximate 
method and is also given in Appendix. One subtle differ- 
ence between the numerical and approximate results is 
zero transmission probability nearby a = ir/2 for exact 
numerical result. However, this physical observation is 
not manifested in the approximate result (dotted line). 
A severe deviation from the dotted line is clearly seen for 
the lines with E > 0.010, which indicates that a excellent 
valley transmission can only be maintained at not very 
large energy. 
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Figure 5: (Color online) The transmission coefficient with 
a fixed momentum q = 0.017/a which is exactly the green 
circle path in Fig. 4. Two lines (blue and red) correspond to 
that of K and K' valleys, respectively. In this case, two line- 
defects are considered and the fixed distance between them is 
N d = 74. 



B. Two line-defects 

The search of the valley polarization material or de- 
vice attracts physics community greatly. Although a pro- 
posal was given recently by scattering off a line- defect pSI 
it looks much difficult for graphene with only one line- 
defect to make electron leave away from the line-defect 
and finally reach the key device, and meanwhile maintain 
at a high valley polarization in this process. Inspired by 
the valley polarization in the line-defect graphene, we ex- 
pect this problem can be resolved by scattering off two 
or several line-defects. 

In Fig. 3, the transmission behavior is mapped as a 
function of the incident angle in the presence of two line- 
defects. Fig. 3(a) shows a lattice model of graphene 
with two line-defects and meanwhile presents a schematic 
interpretation of Nd = 2 which represents the distance 
between two nearest neighbouring line-defects. Using the 
formula (4), the transmission coefficients in the K and K' 
valleys are mapped as a function of the incident angle in 
Fig 3(b) and 3(c). A mirror symmetry with respect to 
a = can be seen clearly between them. Comparing 
to the case of one line-defect, there is almost no distinct 
change for the lines of Nd = 10/30 both in Fig. 3(b) and 
3(c). However, the transmission coefficient starts varying 
for Nd > 70 and in particular some transmission peaks 
can be observed obviously at small incident angles for 
the lines of Nd = 100 and 150. This exciting finding 
makes us fascinated to explore what happens to the two- 
line-defect graphene and why there is a high transmission 
probability at some incident angles. 

In order to demonstrate a whole physical picture, 
we map the transmission coefficient as the function of 
T(q xi q y ) in Fig. 4. The left and right maps, which cor- 
respond to the K and K' valleys, respectively, manifest 
a clear symmetry. Intriguingly, two large transmission 
regions (dark red) can be observed both in the K and K r 




Figure 6: (Color online) Transmission coefficient in only 
K valley is given in the presence of several paralleling line- 
defects: (a) m = 3, (b) m — 5, (c) m = 7, and (d) m — 10. 
Note that here m represents the number of paralleling line- 
defects. In each case, four lines with differen Nd, the dis- 
tance of two nearest neighbouring line-defects, are manifested. 
Here, electron's energy is set as E = 0.01 



valleys. Note that the positive (negative) q x represents 
an incident electron along positive x direction (negative 
x direction). These peaks remind us some quantum in- 
terferences occur in the presence of two line-defects, e.g. 
electron's resonant tunneling phenomena. 

To confirm the origin of the high transmission above, 
we plot two lines selectively in Fig. 5. These two lines 
which belong to K and K' valleys respectively correspond 
to half circle of the green route in Fig. 4, representing 



Q — x/Qx + Qy = 0.017/a. Obviously, the transmission 



peaks at A and A' are from the intersection of the green 
circle route and the high transmission region (dark red). 
When substituting q = 0.017/a and olai A' ~ ±0.057r into 
the formula: A = 1/E = 2/^/3qa and I = Nd/ cos a, we 
find two quantities are equivalent within the physical er- 
ror range. This is to say, the 100% transmission does 
originate from the resonant tunneling, which is often ob- 
served in the system of quantum dot. 33 

According to the quantum theory of wave-particle du- 
ality, the low energy electron in the vicinity of the Dirac 
point is inclined to exhibit undulation characteristics. 
When scattering off two line-defects, a standing wave 
state between two line-defects helps electron resonantly 
tunnel through the barriers of the line-defects. Conse- 
quently, we can observe this perfect transmission phe- 
nomenon in above figures. Using the valley polarization 
formula V = Tk+r*'/ ' ^ e 100% valley polarization can 
not be realized at very small incident angle in the pres- 
ence of two line-defects, e.g. at the positions of the A and 
A' peaks. Nevertheless, an exciting thing is that 100% 
valley polarization can almost be reached when a takes 



Figure 7: (Color online) Contour plots of the transmission 
coefficients in K (left) and K' (right) valleys. In this case, 
six line-defects are considered and the fixed distance between 
arbitrary two neighboring line-defects is Nd = 70. 



the value of [0.27T, 0.37r], a relatively small incident angle. 
At least, the results argue that the high valley transmis- 
sion at a small incident angle is feasible by scattering off 
two line-defects. 



C. Many line-defects 

Inspired by the above results, we expect an excellent 
case could occur to the graphene with several line-defects. 
In Fig. 6, the transmission coefficient only in K valley is 
shown in the presence of several paralleling line-defects: 
(a) m = 3, (b) m = 5, (c) m = 7, and (d) m = 10, where 
m represents the number of paralleling line-defects. In 
Fig. 6(a) with m = 3, almost no obvious change can 
be observed with comparison to the case of m = 2 in 
Fig. 3. Nevertheless, one extra transmission peak still 
appears in the line of Nd = 100 (red). Without question, 
it should be attributed to the effect of multiple interfer- 
ence between line-defects. When setting m = 5 in Fig. 
6(b), except Nd = 30, the transmission coefficients of all 
other lines become active and manifest some peaks at 
some fixed incident angles. For the cases of m = 7 in 
Fig. 6(c) and m = 10 in Fig. 6(d), more transmission 
peaks appears and overlap each other. However, the elec- 
tron transmission coefficient in case of Nd = 30 decreases 
gradually with increasing the number of paralleling line- 
defects. It is because at E = 0.01 electron's wave length 
is considerably larger than the distance of two nearest 
neighbouring line-defects, Nd — 30. 

For sake of concreteness, here we restrict to the case of 
six line-defects and make great efforts to manifest clearly 
a picture of a 100% valley polarization. In Fig. 7, the 
transmission coefficients in two valleys are mapped in the 
panel of (q x ,Qy)- It is clearly shown that comparing to the 
case of two line-defects, more large transmission regions 
(dark red) exist in Fig. 7, the case of six line-defects. 
As explained in Fig. 4 and 5, these high transmission 
regions should likewise be from the resonant tunneling 



Figure 8: (Color online) The transmission coefficient with a 
fixed momentum q = 0.015/a which is exactly the green circle 
path in Fig. 6. Two lines (blue and red) correspond to that 
of K and K valleys, respectively. In this case, six line-defects 
are considered and the fixed distance between arbitrary two 
neighboring line-defects is Nd = 70. 

due to the quantum interference, which is also named as 
Fabry-Prot interference in optics. 

In Fig. 8, the transmission coefficient is plotted as a 
function of the incident angle a, which follows half of the 
green path in Fig. 7. Note that here the electron momen- 
tum q is taken as q = 0.015/a. Obviously, the marking 
A, B, and C transmission peaks should be attributed to 
the three intersections between the green path and the 
large transmission regions (dark red) in Fig. 7. Here, 
a a ~ 0.387T, as ~ 0.157T, and ac ~ — 0.057T illustrate 
that the A transmission peak is attributed to the intrin- 
sic scattering property of a line-defect and the peaks at 
B and C are from the resonant tunneling of the positive 
and negative first standing waves, respectively. 

Much importantly, not only a perfectly resonant trans- 
mission but also a 100% polarization can be obtained 
at the B resonant peak by the polarization formula: 
V = . In a conclusion, we predict that the per- 

fect 100% polarization can be achieved by scattering off 
many line-defects, which origins in the quantum reso- 
nant tunneling phenomenon. Even though only exhibit- 
ing the cases of two and six line-defects here, in principle 
this high transmission and thus the 100% polarization 
can be realized in graphene with arbitrary number of 
line-defects. This feature shed a light on constructing a 
feasible valley polarization resonator by using multiple 
line-defects. 



IV. CONCLUSIONS 

In summary, we study the valley polarization in the 
presence of several paralleling line-defects. Importantly, 
the valley polarization changes with varying the distance 
between the line-defects and the number of the line- 
defects. In particular, the valley polarization can reach to 
the maximum value of one, namely a 100% polarization, 
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at a small incident angle. Further numerical calculations 
manifest that 100% polarization should be attributed to 
the behavior of electron quantum interference, more pre- 
cisely the phenomenon of the standing wave between the 
line-defects. The findings make us get rid of helpless- 
ness in the choose of the electron's incident angle and 
provides a reliable way to realize a controllable valley 
polarization in experiment. With the great advance in 
graphene related technologies, a valley polarization res- 
onator and also a detector of a valley polarization, being 
the key components in valleytronics, are promising to be 
realized in the system of the line-defect graphene. 

In the preparation of this paper, we notice a reference 
on two paralleling line-defect in recentP^and intriguingly, 
nearly a perfect perfect valley polarization is also ob- 
served. This report also illustrates that the valley filter 
of several paralleling line-defects should feasible in exper- 
iment. 
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where Hft, L means a hopping term from the I sub-layer 
in the L supercell to the V sub-layer in the V . Obviously, 

T' njLL rfi njLL 



the relations hold true: Hi 2 L = 



1 2 • 



L 34 
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Thus, the eigen equation [H][C] = E[C] can be written 
as: 
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where E and [C] represents the eigenvalue and eigenstate, 
respectively. For simplicity, we cast the above eigen equa- 
tion into a subspace and consequently obtain the follow- 
ing expression, 



(A.7) 
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APPENDIX 



Here, 



In this appendix, we present a detailed derivation for 
the transmission coefficient T of low energy electron in 
the vicinity of the Dirac point. Note that similar deriva- 
tions can also reference to some previous papers pilMJ 



In the Landauer-Biittiker formula, Eq. Q, the Green's 
function G r can be readily expressed as: G r = 1/(E — 
Hd — — ^r)- Here, Hd is just a two-by-two matrix 
related to T 3 and therefore the key point is to obtain 
the retarded self-energy function, ^ r L ^ R - That can be 
solved by the surface Green's function. In the following, 
we firstly present an analytic method of obtaining the 
surface Green's function G r L and then derive the trans- 
mission coefficient of low energy electron in the vicinity 
of the Dirac point. 

Based on the schematic representation in Fig. 1(c), a 
Hamiltonian of a semi-infinite right lead can be written 



LL 



U L 12 L E- X U^ 
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where T 



E 



nj LL rr — lnj 
H 32 & ft 



23 . Note that other quan- 



tities, such as Wf/s W^ L and Wf 1 L+1 , can be also ob- 
tained by using Dyson equation of Green's function. For 
simplicity in the following derivation, we define some 
symbols as: 



W 



LL 



M/(E 2 - 1), 
= E — ME/(E 2 



1), (A.9) 



= wii L+1 = 1, 



where M takes a form of 
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(A.10) 
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Following the similar derivation processes of the eigen 
equation above, the Green's function can be written read- 
ily as: 



Performing a inverse unitary transformation with re- 
spect to Eq. (A. 13), the surface Green's function of a 



semi-infinite long lead can be expressed as: 
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After a gauge transformation by a unitary matrix, which 
is denoted as: 
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where A = g\ 1 + g\ 2 and B = g\ 1 — g\ 2 . Now, we get an 
exact expression about the surface Green's function of a 
semi-infinite long lead. 

Combining the self energy function and the Green's 
function of the line-defect 

^l/r — r i yi 

G r = l/(E + T 2 f 3 -Xi-Z r R ) (A.17) 



the Ai matrix, all uj$ matrixes, and Green's function in 



Eq. (A.ll) can be transformed into: 
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where mi = 2(1 — cos 0) and m 2 = 2(1 + cos 0). Note that 
the Green's f unctio n must be diagonalized if all coeffi- 
cients in Eq. ( A.ll[ ), namely all uj§, can be diagonalized 
by this unitary transformation. 

As such, Eq. (A.ll) can be decomposed into: 
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where p takes 1 or 2. 

There are many known mathematic method to calcu- 
late an inverse matrix of a tridiagonal matrix. How- 
ever, we refrain from introducing these methods here. 
We use Do to denote the determina nt of the tridiag- 
onal matrix on the left side of Eq. (A. 14) and Di to 



denote the determinant of the tridiagonal submatrix ex- 
cept the former i columns and rows of the tridiagonal 
matrix!^ Then, we can obtain the recursive equations: 
D /D 1 = uj - u 2 D 2 /D 1 and D 1 /D 2 = Zj q - u 2 e D 3 /D 2 . 



In the limit of infinite long lead, corresponding to an 
infinite large tridiagonal matrix in Eq. (A. 14), the re- 
lation D0/D1 = D 2 /Ds maintains. From the formula 
g\ = Di/Dq, we can readily obtain 



dip 



(2; 



2 . 

op 



\2 



■422 



2u5, 



op 

-4-15) 



Substituting the above formula into the Landauer- 
Biittiker formula, Eq. Q, we can express the transmis- 
sion coefficient exactly as: 

T kJ E ) = ^[(Wii)Vil 2 + (Im&1 2 ) 2 |F3| 2 



+8T 2 2 Im^ 1 Im5 1 1 2 cos 2 6'], (A.18) 

where F 1 = 2E- 4rf ^ 2 + 2r 2 cos6», F 2 = (E - At?) 2 - 
Bt 2 (Bt 2 +2t 2 cos0)-t 2 2 , and F 3 = 2£-4t 1 2 ^ 2 -2t 2 cos6». 
In the low energy limit, it is not difficult to find g\ x ~ t, 
E ~ 0, and g\ 2 <C g\i- Thus, the transmission 
coefficient is approximately written asP^ 



012 



TkJE) 



(A.19) 



Note that we still reserve the term of lvag\ 2 m ^1 an d ^2 
because some terms in F\ and F 2 are in the same order 
as lmg[ 2 . 

Assuming incident electron's energy is much smaller 
than the hopping terms t, T\ and r 2 , an approximate 
expression of the transmission coefficient can be written 
as: 



T ky (E) 



t?(1 + cos20)[4E 2 - (1 -mi) 2 ]/2 
2t 1 2 t 2 cos 0(1 — mi) + E 2 + 4t? cos 2 



Following the relations: k x = q x , k y = ±ir/3a + q y , 
E = \/?>qa/2, 6 — k y a « d=7r/3, and tana = q y /q x , 
we adopt the approximation in the low energy limit: 
1 - mi g ± 2Esma, 1 + cos 20 = 1/2, and cos# = 1/2 in 



Eq. (A.20). We can readily obtain 



T K /K'(E) = i(l±sina). 



(A.21) 



op 



That is namely the transmission coefficient of low energy 
electron in the vicinity of two Dirac valleys pHUl 



where uj ep = 1 is adopted. 
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